Enhancement of the Gilbert damping constant due to spin pumping in non-collinear 
ferromagnet / non-magnet / ferromagnet trilayer systems 



o 
o 
o 

as: 
< 

(N 



c 

o 
o 



> 
m 

(N 
m 
cn 

od 
O 

o 



X 



Tomohiro Taniguchi 1 ' 2 , Hiroshi Imamura 2 

1 Institute for Materials Research, Tohoku University, Sendai 980-8577, 
2 Nanotechnology Research Institute, National Institute of Advanced Industrial Science and Technology, 
1-1-1 Umezono, Tsukuba, Ibaraki 305-8568, Japan 
(Dated: February 1, 2008) 

We analyzed the enhancement of the Gilbert damping constant due to spin pumping in non- 
collinear ferromagnet / non-magnet / ferromagnet trilayer systems. We show that the Gilbert 
damping constant depends both on the precession angle of the magnetization of the free layer and 
on the direction of the magntization of the fixed layer. We find the condition to be satisfied to 
realize strong enhancement of the Gilbert damping constant. 
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There is currently great interest in the dynamics of 
magnetic multilayers because of their potential applica- 
tions in non-volatile magnetic random access memory 
(MRAM) and microwave devices. In the field of MRAM, 
much effort has been devoted to decreasing power con- 
sumption through the use of current-induced magnetiza- 
tion reversal (CIMR) 0,H HQ 0, IE 0|- Experimentally, 
CIMR is observed as the current perpendicular to plane- 
type giant magnetoresistivity (CPP-GMR) of a nano pil- 
lar, in which the spin-polarized current injected from the 
fixed layer exerts a torque on the magnetization of the 
free layer. The torque induced by the spin current is 
utilized to generate microwaves. 

The dynamics of the magnetization M in a ferromag- 
net under an effective magnetic field B G ff is described by 
the Landau-Lifshitz-Gilbert (LLG) equation 
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where 7 and ao are the gyromagnetic ratio and the 
Gilbert damping constant intrinsic to the ferromagnet, 
respectively. The Gilbert damping constant is an im- 
portant parameter for spin electronics since the critical 
current density of CIMR is proportional to the Gilbert 
damping constant [8, 9] and fast-switching time magne- 
tization reversal is achieved for a large Gilbert damp- 
ing constant (To| . Several mechanisms intrinsic to ferro- 
magnetic materials, such as phonon drag 11] and spin- 
orbit coupling [l2| . have been proposed to account for 
the origin of the Gilbert damping constant. In addition 
to these intrinsic mechanisms, Mizukami et a/. [HI, Il4j 
and Tserkovnyak et al. [TH, EE showed that the Gilbert 
damping constant in a non-magnet (N) / ferromagnet (F) 
/ non-magnet (N) trilayer system is enhanced due to spin 
pumping. Tserkovnyak et aL[l7| also studied spin pump- 
ing in a collinear F/N/F trilayer system and showed that 
enhancement of the Gilbert damping constant depends 
on the precession angle of the magnetization of the free 
layer. 

On the other hand, several groups who studied CIMR 
in a non-collinear F/N/F trilayer system in which the 
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FIG. 1: (Color online) The F/N/F trilayer system is schemat- 
ically shown. The magnetization of the Fi layer (mi) pro- 
cesses around the z-axis with angle 6 and angular velocity ui. 
The magnetization of the F2 layer (ma) is fixed with tilted 
angle p. The precession of the magnetization in the Fi layer 
pumps spin current JP um P into the N and F2 layer, and creates 
the spin accumulation /in in the N layer. The spin accumu- 



lation induces the backflow spin current I, 



back(i) 



(< = 1,2). 



magnetization of the free layer is aligned to be perpen- 
dicular to that of the fixed layer have reported the reduc- 
tion of the critical current density [E IE Q- Therefore, it 
is intriguing to ask how the Gilbert damping constant is 
affected by spin pumping in non-collinear F/N/F trilayer 
systems. 

In this paper, we analyze the enhancement of the 
Gilbert damping constant due to spin pumping in non- 
collincar F/N/F trilayer systems such as that shown in 
Fig. [TJ Following Refs. [IE HE HE El , we calculate the 
spin current induced by the precession of the magnetiza- 
tion of the free layer and the enhancement of the Gilbert 
damping constant. We show that the Gilbert damping 
constant depends not only on the precession angle 9 of 
the magnetization of a free layer but also on the angle p 
between the magnetizations of the fixed layer and the pre- 
cession axis. The Gilbert damping constant is strongly 
enhanced if angles and p satisfy the condition 6 = p or 
= tt — p. 

The system we consider is schematically shown in Fig. 
[JJ A non-magnetic layer is sandwiched between two fer- 
romagnetic layers, Fi and F2. We introduce the unit 
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vector mi to represent the direction of the magnetiza- 
tion of the i-th ferromagnetic layer. The equilibrium 
direction of the magnetization rti\ of the left free fer- 
romagnetic layer Fi is taken to exist along the z-axis. 
When an oscillating magnetic field is applied, the magne- 
tization of the Fi layer precesses around the z-axis with 
angle 9. The precession of the vector mi is expressed 
as mi = (sin 9 cos uit, sin 9 sin u>t, cos 6), where ui is the 
angular velocity of the magnetization. The direction of 
the magnetization of the F 2 layer, m 2 , is assumed to be 
fixed and the angle between m 2 and the z-axis is repre- 
sented by p. The collinear alignment discussed in Ref. 
[ijj corresponds to the case of p = 0, n. 

Before studying spin pumping in non-collinear sys- 
tems, we shall give a brief review of the theor y o f 
spin pumping in a collinear F/N/F trilayer system [17| . 
Spin pumping is the inverse process of CIMR where the 
spin current induces the precession of the magnetization. 
Contrary to CIMR, spin pumping is the generation of 
the spin current induced by the precession of the mag- 
netization. The spin current due to the precession of the 
magnetization in the Fi layer is given by 
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where g^l [ s a mixing conductance [H, Gil and fi is the 
Dirac constant. Spins are pumped from the Fi layer 
into the N layer and the spin accumulation /xn is cre- 
ated in the N layer. Spins also accumulate in the Fi 
and F 2 layers. In the ferromagnetic layers the trans- 
verse component of the spin accumulation is assumed to 
be absorbed within the spin coherence length defined as 



A, 



fcp |, where kp is the spin-dependent 
Fermi wave number of the i-th ferromagnet. For fer- 
romagnetic metals such as Fe, Co and Ni, the spin co- 
herence length is a few angstroms [20| . Hence, the spin 
accumulation in the i-th ferromagnetic layer is aligned to 
be parallel to the magnetization, i.e., /xf ; — pp i mi. The 
longitudinal component of the spin accumulation decays 
on the scale of spin diffusion length, A^j , which is of the 
order of 10 nm for typical ferromagnetic metals [2l| . 

The difference in the spin accumulation of ferromag- 
netic and non-magnetic layers, A/Xj = /xn — MFi m i(* = 
1,2), induces a backflow spin current, j^ ack ( 2 ^ flowing 
into both the Fi and F 2 layers. The backflow spin cur- 



rent I, 



back(i) 



is obtained using circuit theory [181 ] as 
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where g'' and g^ are the spin-up and spin-down con- 
ductances, respectively. The total spin current flowing 
out of the Fi layer is given by I c * ch = J pump - J 6 back W 
[l7j |. The spin accumulation /ip i in the Fj layer is ob- 
tained by solving the diffusion equation. We assume 
that spin-flip scattering in the N layer is so weak that 



we can neglect the spatial variation of the spin current 



within the N layer, I\ 



cxch 



r back(2) 



The torque T\ 



acting on the magnetization of the Fi layer is given by 



Ti = I° xch - (mi • I™>i = mi x (J| xch 
the collinear system, we have 



x mi). For 
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where v — (gM-l _ g*)/(g^ + g*) is the dimensionless 
parameter introduced in Ref. [171 ] . The Gilbert damping 
constant in the LLG equation is enhanced due to the 
torque Ti as «o ~~ * a o + ol with 
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where <7l is the Lande g-factor, /ie is the Bohr magneton, 
dpi is the thickness of the Fi layer and S is the cross- 
section of the Fi layer. 

Next, we move on to the non-collinear F/N/F trilayer 
system with p = 7r/2, in which the magnetization of the 
F 2 layer is aligned to be perpendicular to the z-axis. Fol- 
lowing a similar procedure, the LLG equation for the 
magnetization M\ in the Fi layer is expressed as 

dMi . . „ , 7off, . u Mi dM x 

-^ r = - 7cff Mixi3 cfr+ — (ao+aj — x—, (6) 

where 7 c g and a are the effective gyromagnetic ratio 
and the enhancement of the Gilbert damping constant, 
respectively. The effective gyromagnetic ratio is given by 

ft gi^^g^v cot9 cosip sinuit\ 1 . . 
7cff= H 8nMd Fl Se ) ' (?) 

where cosV' = sin 9 cos = m\- m 2 and 

e = 1 — v 2 cos 2 ip — ^(cot 2 9 cos 2 ip — sin 2 ip + sin 2 u>t) . (8) 

The enhancement of the Gilbert damping constant is ex- 
pressed as 
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v cot 2 9 cos 2 ip 
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It should be noted that, for non-collinear systems, both 
the gyromagnetic ratio and the Gilbert damping constant 
are modified by spin pumping, contrary to what occurs 
in collinear systems. The modification of the gyromag- 
netic ratio and the Gilbert damping constant due to spin 
pumping can be explained by considering the pumping 
spin current and the backflow spin current [See Figs. [2]^ a) 
and[^b)]. The direction of the magnetic moment car- 
ried by the pumping spin current /P um P is parallel to 
the torque of the Gilbert damping for both collinear and 
non-collinear systems. The Gilbert damping constant is 
enhanced by the pumping spin current JP ump . On the 
other hand, the direction of the magnetic moment carried 



3 



l s pump 



F1 | s back(1) 



P ^ |M2 

,back(1) F2 



(b) 



backflow 



precession 




Mi 

Gilbert damping 

FIG. 2: (Color online) (a) Top view of Fig. □ The dotted 
circle in Fi represents the precession of magnetization Mi 
and the arrow pointing to the center of this circle represents 
the torque of the Gilbert damping. The arrows in JP um P and 
jback(i) re p regen t ^he magnetic moment of spin currents, (b) 
The back flow j^ 3 ^ 1 ) ij as components aligned with the di- 
rection of the precession and the Gilbert damping. 



by the backflow spin current j^^M depends on the di- 
rection of the magnetization of the F2 layer. As shown in 
Eq. © , the backflow spin current in the F2 layer jback(2) 
has a projection on m^. Since we assume that the spin 
current is constant within the N layer, the backflow spin 
current in the Fi layer j^^C 1 ) a i so h as a projection on 
m,2. For the collinear system, both JP um P and jkack(i) 
are perpendicular to the precession torque because m<i 
is parallel to the precession axis. However, for the non- 



collinear system, the vector I, 



back(l) 



has a projection on 
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dt 



MA 



(10) 



The estimated value of a is of the order of 0.001. How- 
ever, we cannot neglect a since it is of the same order 
as the intrinsic Gilbert damping constant ao [22I [23[ . 

Experimentally, the Gilbert damping constant is mea- 
sured as the width of the ferromagnetic resonance (FMR) 
absorption spectrum. Let us assume that the Fi layer has 
no anisotropy and that an external field B C7< ± = BqZ is 



applied along the z-axis. We also assume that the small- 
angle precession of the magnetization around the z-axis 
is excited by the oscillating magnetic field B\ applied in 
the sy-plane. The FMR absorption spectrum is obtained 
as follows [24j : 



dt 



ajMfl 2 B 2 



(7S0 - fi) 2 + ( ai B y 



(11) 



the precession torque, as shown in Fig. [DJb). Therefore, 
the angular momentum injected by jb ack W modifies the 
gyromagnetic ratio as well as the Gilbert damping in the 
non-collinear system. 

Let us estimate the effective gyromagnetic ratio using 
realistic parameters. According to Ref. [l7j . the con- 
ductances g'* and g* for a Py/Cu interface are given 
by g^/S = 15[nm~ 2 ] and v ~ 0.33, respectively. The 
Lande g- factor is taken to be <?l = 2.1, magnetization is 
AttM = 8000 [Oe] and thickness d Fl = 5[nm]. Substitut- 
ing these parameters into Eqs. |J7J| and {SJ, one can see 
that |7cff/7— 1| — 0.001. Therefore, the LLG equation 
can be rewritten as 



where fi is the angular velocity of the oscillating mag- 
netic field, T — 2tt/Q and a — ctQ + a . Since a is very 
small, the absorption spectrum can be approximately ex- 
pressed as P oc a?o + (a ) and the highest point of the peak 
proportional to (l/(ao + a )), where (a') represents the 
time-averaged value of the enhancement of the Gilbert 
damping constant. In Fig. El^a), the time-averaged value 
(a ) for a non-collinear system in which p — tt/2 is plot- 
ted by the solid line as a function of the precession an- 
gle 0. The dotted line represents the enhancement of 
the Gilbert damping constant a for the collinear system 
given by Eq. ([5]). The time- averaged value of the en- 
hancement of the Gilbert damping constant (a') takes 
its maximum value at 9 — 0, 7r for the collinear system 
(p = 0, 7r). Contrary to the collinear system, (a'} of the 
non-collinear system in which p = tt/2 takes its maxi- 
mum value at 9 — tt/2. 

As shown in Fig. HJb), the backflow spin current 
gives a negative contribution to the enhancement of the 
Gilbert damping constant. This contribution is given by 
the projection of the vector jb ack W onto the direction 
of the torque of the Gilbert damping, which is repre- 
sented by the vector mi x rhi. Therefore, the condition 
to realize the maximum value of the enhancement of the 
Gilbert damping is satisfied if the projection of j^ ack ^ 
onto mi x rhi takes the minimum value; i.e., 6 = p or 

= TT — p. 

We can extend the above analysis to the non-collinear 
system with an arbitrary value of p. After performing the 
appropriate algebra, one can easily show that the LLG 
equation for the magnetization of the Fi layer is given 
by Eq. © with 



7eff = 7 



3lMbS ^sinpsinwi(cot#cos'0 — csc#cosp) 
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(13) 



where cos?/> = sin 9 sin p cos uit + cos 9 cos p = mi ■ m-2 
and 

e =1 — v 2 cos 2 4> 

— ^{(cot 9 cos ip — esc 9 cos p) 2 — sin 2 ip + sin 2 p sin 2 ut} . 

(14) 

Substituting the realistic parameters into Eqs. (fT2|) and 
([H]), we can show that the effective gyromagnetic ratio 
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FIG. 3: (Color online) (a) The time-averaged value of the en- 
hancement of the Gilbert damping constant a is plotted as a 
function of the precession angle 8. The solid line corresponds 
to the collinear system derived from Eq. ((9j). The dashed 
line corresponds to the non-collinear system derived from Eq. 
(0). (b) The time-averaged value of the enhancement of the 
Gilbert damping constant a of the non-collinear system is 
plotted as a function of the precession angle and the an- 
gle p between the magnetizations of the fixed layer and the 
precession axis. 



7etr can be replaced by 7 in Eq. (J6j) and that the LLG 
equation reduces to Eq. (fTT))) . Figure [31(b) shows the 
time-averaged value of the enhancement of the Gilbert 
damping constant (a ) of Eq. (flU|) . Again, the Gilbert 
damping constant is strongly enhanced if angles 9 and p 
satisfy the condition that 6 = p or 9 = tt — p. 



In summary, we have examined the effect of spin 
pumping on the dynamics of the magnetization of mag- 
netic multilayers and calculated the enhancement of the 
Gilbert damping constant of non-collinear F/N/F trilayer 
systems due to spin pumping. The enhancement of the 
Gilbert damping constant depends not only on the pre- 
cession angle 9 of the magnetization of a free layer but 
also on the angle p between the magnetizations of the 
fixed layer and the precession axis, as shown in Fig. [3Jb) . 
We have shown that the 9- and p-dependence of the en- 
hancement of the Gilbert damping constant can be ex- 
plained by analyzing the backflow spin current. The con- 
dition to be satisfied to realize strong enhancement of the 
Gilbert damping constant is 9 = p or 9 = n — p. 
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